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Time-dependent nonequilibrium Green's functions are used to study electron trans- 
port properties in a device consisting of two linear chain leads and a time-dependent 
interleads coupling that is switched on non-adiabatically. We derive a numerically ex- 
act expression for the particle current and examine its characteristics as it evolves in 
time from the transient regime to the long-time steady-state regime. We find that just 
after switch-on the current initially overshoots the expected long-time steady-state 
value, oscillates and decays as a power law, and eventually settles to a steady-state 
value consistent with the value calculated using the Landauer formula. The power- 
law parameters depend on the values of the applied bias voltage, the strength of 
the couplings, and the speed of the switch-on. In particular, the oscillating tran- 
sient current decays away longer for lower bias voltages. Furthermore, the power-law 
decay nature of the current suggests an equivalent series resistor-inductor-capacitor 
circuit wherein all of the components have time-dependent properties. Such dynam- 
ical resistive, inductive, and capacitive infiuences are generic in nano-circuites where 
dynamical switches are incorporated. We also examine the characteristics of the 
dynamical current in a nano-oscillator modeled by introducing a sinusoidally modu- 
lated interleads coupling between the two leads. We find that the current does not 
strictly follow the sinusoidal form of the coupling. In particular, the maximum cur- 
rent does not occur during times when the leads are exactly aligned. Instead, the 
times when the maximum current occurs depend on the values of the bias potential, 
nearest-neighbor coupling, and the interleads coupling. 
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^ ^eduar do . cuansingOgmail . com| ^- ^elelgOnus . edu . s"g| 



1 



I. INTRODUCTION 

The further miniaturization of electronic 
devices wiU eventuaUy lead to molecular 
electronics wherein particles pass through 
molecular-scale devices whose constituent 
molecules may have been manipulated and 
synthetically assembled or created.- Molec- 
ular transistors, in particular, are of signifi- 
cant practical interests and whose successful 
implementations are currently actively being 
pursued. Several theoretical models of the 
transistor have been proposed^^^ and exper- 
imental successes have also been reported.^ 
A related molecular-scale device, the molec- 
ular switch, has also garnered significant in- 
terests because of the switch's important role 
in circuit design and architecture. A theo- 
retical model of the switch makes use of a 
mechanism that involves the reversible dis- 
placement of an atom in a molecular wire 
through the application of a gate voltage.^ 
In addition, experimental realizations of the 
atomic switch include mechanisms involving 
the reversible transfer of an atom between 
two leads,- the dynamical onset of single- 
atom contact between leads'" and the manipu- 
lation of atomic bonds using a dynamic force 
microscope.^ Having a dynamical switch in 
a nano-circuit, however, necessitates the ap- 
pearance of time-dependent behavior, partic- 



ularly during the transient regime just after 
switch-on. The circuit transitions from being 
disconnected into connected in a short time 
and the current does not instantly switches 
into the steady-state value upon connection. 
It is therefore informative to know the char- 
acteristics of the current just after switch- 
on and during the transient regime, and 
determine how this current approaches the 
steady-state value. In this work, we intro- 
duce a model device representing a system 
wherein the current can be dynamically tog- 
gled on and off. There are several theoret- 
ical approaches in treating time-dependent 
quantum transport. Among these include 
time-dependent density functional theory,- 
propagating the Kadanoff-Baym equations,— 
Floquet theory,^^ path-integral techniques,— 
and the density matrix renormalization group 
method.—"— In this work we choose to 
use time-dependent nonequilibrium Green's 
functions (TD-NEGF) because, as is shown 
in this paper, dynamically toggling the cou- 
pling between the leads is straightforward us- 
ing the technique, even without the assump- 
tion of weak coupling, and the resulting ex- 
pression for the time-dependent current is nu- 
merically exact. 

The device we examine consists of a semi- 
infinite linear chain of atoms, i.e., the left 
lead, which is stationary and another semi- 
infinite linear chain of atoms, i.e., the right 
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FIG. 1. (Color online) (a) An illustration of the 
nano-relay and nano-oscillator. The left lead is 
firmly on a stationary substrate and the right 
lead is on a rotatable disk. The left and right 
leads align when the disk is rotated clockwise 
to the dash line, (b) The device can be repre- 
sented by two semi-infinite leads connected by 
a time- varying coupling v^^{t) switched on at 
t = 0. In both leads, the on-site energy is e and 
the nearest-neighbor hopping parameter is v. (c) 
An equivalent resistor-inductor-capacitor circuit 
with dynamical properties. 

lead, that is on a rotatable disk. An illus- 
tration of the device is shown in Fig. [l](a). 
When the disk is rotated clockwise to the 
dashed line, the left and right leads align and 
conduction occurs. We model how the two 
leads couple in two ways: as an abrupt Heav- 
iside step function and as a hyperbolic tan- 
gent that gradually progresses in time. For 
such cases, the off state has no current flow- 
ing across the leads. In addition, we can also 
model an oscillator by swinging the disk back 
and forth across the dashed line. The ro- 



tatable disk can also be replaced by a gate 
voltage, located on top of the right lead, that 
could rotate the right lead to its desired posi- 
tion. This latter configuration has previously 
been examined and, in particular, the steady- 
state transport properties of its on and off 
states have been studied^. In this paper, 
however, we study the full time-dependent 
transport properties of the device. Numeri- 
cally exact expressions for the current and the 
needed nonequilibrium Green's functions are 
derived. We then show that just after switch- 
on and during the transient regime the cur- 
rent initially overshoots the expected value of 
the steady-state current and then oscillates 
around the steady-state value while decay- 
ing as a power law. Such a power-law decay- 
ing transient current suggests the appearance 
of dynamical resistance, inductance, and ca- 
pacitance in the system during the transient 
regime. Power-law decaying currents have 
also recently been predicted to appear in a 
system containing a quantum dot channel^ 
and in the anisotropic Kondo model.— 

II. MODEL AND METHOD 

We implement time-dependent nonequi- 
librium Green's functions techniques to in- 
vestigate the dynamical transport proper- 
ties of particles traversing through a device 
with time- varying components. To deter- 



mine the dynamical transport properties, a 
TD-NEGF approach can be used that uti- 
hzes either two-time^^ii^ Green's functions, or 
double-energy-^^ Green's functions, or Green's 
functions that depend on one time and one 
energy variables.— In a transistor with 
source-channel-drain and top gate configura- 
tion, for example, the device can be mod- 
eled by a Hamiltonian constructed using den- 
sity functional theory2ii2^ or tight-binding 
theorySi and the dynamical transport prop- 
erties of small channels are calculated using 
TD-NEGF. For devices with larger channels, 
a self-consistent calculation based on the 
Poisson equation and TD-NEGF can be done 
to determine the consistent dynamical poten- 
tial and charge density in the channel. In 
this work, in contrast, we study a device con- 
sisting only of two leads, i.e., there is no chan- 
nel between the leads. The time-dependence 
comes from non-adiabatically toggling on the 
coupling between the leads. We use TD- 
NEGF to derive a numerically exact expres- 
sion for the dynamical current and investi- 
gate how the devices we call a nano-relay 
and a nano-oscillator respond to time- varying 
influences. This approach has recently also 
been used in the study of dynamical heat 
transport in a thermal switch^^. 

We model the system by the total Hamil- 
tonian H = + + H^^, where is 
for the left lead on the stationary substrate. 



is for the right lead on the rotatable 
disk, and H^^ includes the dynamic inter- 
leads coupling. In the leads, particles follow 
the tight-binding Hamiltonian 

k kj 

(1) 

where c^^ and are particle creation and an- 
nihilation operators at the kth site in the a 
lead. Ck is the on-site energy at site k and Vkj 
is the hopping parameter between nearest- 
neighbor sites k and j (see Fig. [11(b)). The 
sums are over all sites in the a lead. H^^ in- 
cludes the time-dependent component of the 
total Hamiltonian and is of the form 

^'^ = E(<^'^?+4's^'^^)' (2) 

kj 

where v^^{t) = v^j"{t) is the time-dependent 
coupling between the left and right leads and 
is switched-on at t = 0. Only the right-most 
site of the left lead, i.e., the site labeled in 
Fig. [11(b), can couple to the left-most site of 
the right lead, i.e., the site labeled 1 in the 
flgure. 

The current can be determined by not- 
ing how the number operator, A^^ = 
^^c^^c^, changes with time, i.e., I^{t) = 
—q l^dN^ /d£)^ where q is the electron charge. 
Deflning the two-time lesser Green's function 
as 

G'J"'<(ti,t2) = ^(4^(t2)c»), (3) 



we can write the electron current flowing out of the right lead as 
I 
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Similar steps can be done to determine the steady-state regimes, we employ a contour 
current flowing out of the left lead. We flnd that begins at t = when the interleads cou- 
I^{t) = —I^{t) and therefore, current is al- pling has just been switched on, then pro- 
ways conserved at each instant of time t. ceeds to time t where we want to determine 
We deflne the contour-ordered Green's the current, and then goes back to time t = 0. 
function the interaction picture the contour-ordered 

Green's function can be expanded in powers 

Gft{ri,r2) = -^(t, cf{Ti)c];\T2)y (5) of i/h. Applying Wick's theorem to the re- 
sulting expansion and then using Langreth's 

where Tc is the contour-ordering opera- theorem and analytic continuation,— we get a 

tor and ti and T2 are variables along the numerically exact expression that includes all 

contour.^'' Since we want to calculate the terms in the expansion for the lesser Green's 

current for both the steady-state and non- function: 
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(6) 

where the advanced and retarded Green's functions are 

Gf'Hh, h) = Gffih, h) + E r Gf^i(ti,t')vl:^nit')G^iHt', h), C = r, a. (7) 

mn 

The first-order retarded and advanced Green's functions are 

Gffiti, h) = f dt' gff{t, - t')vf^{t')g'^i{t' - h), (8) 
Jo 



and the first-order lesser Green's function is 

G^fituh) = 1^ dt' [gffih - t')vf^{t')g]:i:'{t' - h) + gfj<{h - t')vf{t')g]:':{t' -h)}, 

(9) 

I 

where the ^-^(t), and g<{t) are the III. RESULTS AND DISCUSSION 

retarded, advanced, and lesser free-leads 

Green's functions, respectively. Time- Firstly, we examine the impact of the 

translation invariance is satisfied by the switch-on speed to current characteristics, 

free leads and therefore, their correspond- The functional form of the interleads cou- 

ing Green's functions can be calculated in pling, v^'^it), describes how the device is 

the energy domain using the techniques switched on. We examine two types of 

of steady-state NEGF.-^^ The integrals in switch-ons: an abrupt Heaviside step func- 

Eqs. (!HD and © are then determined us- tion switch-on and a gradually progressing 

ing the extended Simpson's rule algorithm hyperbohc tangent switch-on of the form 

for numerical integration.^ Furthermore, the ^^^^(^) = v^^tanh{udt), where cod is the driv- 

expressions for the advanced and retarded ing frequency. The step function switch-on is 

Green's functions in Eq. (CD are in the actually the hmit when ^ oo of the hy- 

form of a Fredholm equation of the sec- perbohc tangent switch-on. Furthermore, we 

ond kind and can be solved by discretizing on-site energy e = 0. The left and 

the time integrals, based on the extended right leads have temperature = = 

Simpson's rule algorithm, and performing a 300 K and chemical potential /il = and 

matrix inversion using LU (lower-triangular /^R = - qVb, where we set the Fermi en- 

and upper-triangular) decomposition.-^! The ergy = 0. The bias potential H is apphed 

lesser Green's function can then be calculated to the right lead and the bias potential energy 

from the retarded and advanced Green's i^ written as Ut, = qVb. 

function by applying the extended Simpson's Figure M shows the current flowing out of 

rule algorithm to numerically solve the inte- the left lead for the step function and hy- 

grals in Eq. dHj). perbolic tangent switch-on with driving fre- 
quencies u!d = 0.25 [1/t] and 0.5 [1/t], where 
[1/t] = 10-^^ rad/s. The steady-state val- 
ues of the current calculated separately us- 
ing the Landauer formula for a linear chain 




time [fs] 



FIG. 2. (Color online) The current clS cl function 
of time when the interleads couphng v^^{t) is 
switched on as a step function at t = (black 
hues), graduahy as a hyperbohc tangent with 
driving frequency cjd = 0.5 [1/t] (red hnes), and 
(jjd = 0.25 [1/t] (green hnes), where [1/t] = 
10-*^^ rad/s. The bias potentials are (a) Ui, = 
0.5 eV, (b) Ut = 0.3 eV, and (c) Ut = 0.1 eV. 
The (blue) dashed lines are the values of the 
steady-state current. The hopping parameter in 
the leads is v = —2.7 eV. 

with unit transmission^^ are also shown in 
Fig. [2] as dashed lines. During the times 
just after the switch-on, the current rapidly 
increases and overshoots the expected long- 
time steady-state value. It then oscillates 
and decays in time, eventually settling to the 
steady-state value. It can be seen that as the 
driving frequency is decreased the amplitude 
of the oscillations also decreases. In addition, 
the peaks are displaced to later times because 
of the more gradual progress of the interleads 
coupling v^^{t). In Fig. [2] we also see the de- 



pendence of the decay time of the transient 
current to the applied bias potential and the 
speed of the switch-on. The higher bias re- 
sults in a faster decay time for the oscillating 
transient current. 

The transient current oscillates and de- 
cays in time until it settles to a steady- 
state value. During the transient regime 
the strength of the interleads coupling dy- 
namically changes resulting in particles tem- 
porarily accumulating at the left and right 
sides of that coupling. Although we do not 
explicitly consider Coulomb interactions be- 
tween charges, the temporary accumulation 
of charges at the sides of the interleads cou- 
pling, together with the distance between the 
accumulated charges and the applied bias 
voltage across the interleads coupling, can 
be regarded to generate a temporary dy- 
namical capacitance. Similarly, a dynami- 
cal inductance may arise because the tran- 
sient current is varying in time. By consider- 
ing possible equivalent circuit combinations 
and performing least-squares fitting to the 
envelope of the decaying transient current we 
find that the decay closely follows a power 
law, indicating an equivalent series resistor- 
inductor-capacitor (RLC) circuit whose com- 
ponents have time-dependent properties. In 
a previous study using semiclassical Boltz- 
mann transport theory on quantum wires, it 
is found that the wire can be modeled by an 



equivalent series RLC circuit. A series RLC itance vary in time, the resulting transient 
circuit consisting of components with con- current can oscillate and could decay as a 
stant resistance, inductance, and capacitance power law. Making therefore such an anal- 
results in a transient current whose envelope ogy to the quantum device we are examining 
decays as an exponential function. However, (see Fig. [l](c)), applying Kirchhoff's law to a 
when the resistance, inductance, and capac- series RLC circuit with time-varying compo- 
nents leads to the equation 



(fl fR ldL\dI fldR ^\t_ Q dC 



r 



where / is the time-dependent current and ( IT2l ) we get the coupled equations 
through the circuit, R = R{t) is the resis- ^ I dL 2a 

tance, L = L{t) is the inductance, C = C{t) L ~^ L dt t ' 

is the capacitance, and Q = I dt is the — 7^ ^^t^ (^2) 



time-dependent charge accumulating at the Q dC 2 



L dt LC 



capacitor. Furthermore, the power law fits 
imply that the current is of the form which can be solved to determine how 

L(t), and C{t) vary in time for specific values 

/(t) = /or"sinM + (/)) + /o, (11) of "and a;. 

The power-law exponent a determines 

how fast the transient current decays until it 

where a is the power-law exponent de- 
reaches the steady-state value. In Fig. [3](a), 

termined from the fits, u is the time- 
it can be observed that by increasing the 

independent frequency of oscillation of the 

bias potential the value of a also increases, 

transient current, (j) is the phase determined 

thereby speeding up the decay of the tran- 

from initial conditions, and Iq is the time- 

sient current, a and the power-law coeffi- 

independent steady-state current. Taking the 

cient Iq actually also follow power-law fits 

time derivative of Eq. f fTTI) twice, we find 

when the bias potential is varied, as can be 
^ seen in Fig. [31 It suggests a = ao uP and 

d^I 2a(i/ 2r_ 2r /-,9\ 

~d^ ~^ ~r~dt ~^ ~ 0' V ) Iq = Iqo , where ao and Iqq are indepen- 
dent of [4. The power-law exponents /3 and 
where cj^^ = ^2_^ ^(^-1) ^ Comparing Eqs. f fTOl ) 7 determine how fast a and /q, respectively, 
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FIG. 3. (Color online) (a) The power-law ex- 
ponent a and (b) coefficient Iq as functions of 
the bias potential. The (green) dots are for the 
Heaviside step function interleads coupling. The 
(red) triangles are for the hyperbolic tangent in- 
terleads coupling with ujd = 0.5 [1/t]. The (blue) 
squares are for cjd = 0.25 [1/t]. The dashed 
lines are power-law fits to the corresponding 
data points. The amplitude of the couplings are 
y — ^LR — —2.7 eV. Error bars arise from the 
difference between values from the top and bot- 
tom envelopes. 

change when the bias potential is varied. In 
Tabled] we show how the values of the power- 
law fitting parameters change when [/^ is var- 
ied. 

From Table [U we see that the values of 
Iqo and 7 are independent of the speed of 
the switch-on. In addition, the exponent 7 is 
about one. These imply that Iq increases lin- 
early with the bias potential and is consistent 
with the identification that Iq is the steady- 
state current. For the exponent a, we find 
that as the speed of the switch-on is increased 



TABLE I. Values from the power-law fits to a 
and Iq when the interleads coupling is in the 
form of a hyperbolic tangent and a step func- 
tion. The couplings are v = —2.7 eV and 
v^^ = —2.7 eV. The dimension of ao is (1/eV)^ 
and /qo is /xA/(eV)^. 
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0.392 
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0.307 


0.332 


40.214 


0.983 



the coefficient ao also increases but the ex- 
ponent /3 decreases. The increasing ao sug- 
gests that for a given bias potential, the faster 
switch-on results in a faster decay of the tran- 
sient current. Since the slightly decreasing /3 
is still positive, increasing the bias potential 
still speeds up the decay of the transient cur- 
rent. As a result, when the device is operated 
under low bias its current suffers oscillations 
and overshootings longer than when it is op- 
erated under higher bias. Furthermore, it can 
be observed that the power-law parameters 
ao and /3 control the speed of decay of the 
transient current. The values of these two 
parameters vary depending on the speed of 
the switch-on. If we want the system to have 
a fast decaying transient current, then our re- 
sults indicate that we need a switch-on that 
is as fast as possible. 

Next, we investigate the effects of vary- 
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FIG. 4. (Color online) The current function 
of time for a step function switch-on. (a) The 
couplings have values v = v^^ = —4.0 eV (red 
line), V = v^^ = —2.7 eV (black line), and v = 
v^^ = -2.0 eV (green line). The (blue) dashed 
line is the value of the steady-state current, (b) 
The hopping parameter is fixed at v = —2.7 eV 
and the interleads coupling is varied with values 
i;LR ^ _2.7 eV (black line), v^^ = -2.4 eV 
(violet line), and = —2.1 eV (orange line). 
The bias potential for both plots is Ub = 0.3 eV. 

ing the values of the hopping parameter v 
and the interleads coupling v^^. The val- 
ues of these tight-binding parameters depend 
on the material used and varying them effec- 
tively means that we change the material we 
use for the device. The parameters can be 
varied separately or they can be varied while 
maintaining v = v^^. Firstly, we consider the 
latter case. In order to determine only the 
effects of varying the couplings, and not the 
effects of the speed of the switch-on, we em- 
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ploy the step function switch-on. The current 
as a function of time is shown in Fig. [4](a). 
Since v = v^^^ the long-time steady-state val- 
ues of the current can be calculated from the 
Landauer formula with a transmission coef- 
ficient T = 1 (no scattering involved). This 
steady-state value is shown in Fig. [4](a) as 
a dashed line. We examine coupling values 
V = v^^ = -2.0 eV, -2.7 eV, and -4.0 eV. 
The bias potential is fixed at C/5 = 0.3 eV. 
We find that as the couplings become more 
negative the frequency of oscillation of the 
transient current increases. The decrease in 
the value of the couplings imply that the en- 
ergy needed for the particle to hop from one 
site to a neighboring site is decreased. This 
frees up the particle, thereby allowing higher 
oscillation frequencies in the transient cur- 
rent. For long times after the transient os- 
cillations have decayed away, the steady-state 
value of the current is independent of the spe- 
cific value of the couplings. 

Furthermore, we examine the effects of 
varying v and v^^ separately. When v^^ is 
different from the interleads distance is dif- 
ferent from the nearest-neighbor distance be- 
tween sites in the leads. This results in a 
potential barrier that is different at the in- 
terleads coupling and thus, a particle moving 
from the left lead scatters at the interleads 
coupling. Figure [4](b) shows the current as 
a function of time when the nearest-neighbor 



hopping parameter is set at 'U = — 2.7 eV and 
we vary v^^ to values —2.7 eV, —2.4 eV, and 
—2.1 eV. We do not consider > \v\ 

values because that would imply a shorter 
interleads distance than the natural nearest- 
neighbor distance, represented by 'U, in the 
leads. In contrast, decreasing |'^^^| increases 
the potential barrier at the interleads cou- 
pling, and thus implying a longer interleads 
distance, and results in the reduction in the 
amplitude of the oscillating transient cur- 
rent. From Fig. [11(b), we also see that the 
peaks are slightly shifted to later times. In 
addition, the long-time steady-state current 
slightly decreases when \v^^\ is decreased. If 
we want to calculate the steady-state cur- 
rent using the Landauer formula, we would 
find that the transmission coefficient is re- 
duced when v^^ is different from v because 
of the scattering occuring at the interleads 
coupling. 

The oscillating and decaying transient cur- 
rent when V and v^^ are varied also follow a 
power law. In Figs. [5](a) and [5](b) the val- 
ues of V and v^^ are varied together while 
in Figs. [5](c) and [5](d) the values are varied 
separately. Figure [5](b) shows that the val- 
ues of Iq is the same for the cases examined. 
Compared to Fig. [5](d), we see that the val- 
ues of Iq are slightly different. Identifying Iq 
as the steady-state current, we thus confirm 
that it is the same whenever v = v^^. On the 
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FIG. 5. (Color online) The exponent a and co- 
efficient Iq when v and v^^ are varied. In (a) 
and (h)^ V — v^^ = —2.0 eV for the (red trian- 
gles), V = v^^ = —2.7 eV for the (black dots), 
and V = v^^ = —4.0 eV for (blue squares). In 
(c) and (d), we set = — 2.7 eV and vary v^^ 
to values —2.1 eV (green diamonds), —2.4 eV 
(orange inverted triangles), and —2.7 eV (black 
dots). The dashed lines are the power-law fits to 
the corresponding data points. 

other hand, the scattering that happens at 
the interleads coupling when v^^ is different 
from V affects the value of the steady-state 
current. Moreover, the plots of the exponent 
a and coefficient Iq as functions of the bias 
potential can also be fitted to power laws. 
As shown in Table [TTl when v = v^^^ de- 
creasing the value of the couplings decreases 
both the coefficient ao and the exponent /3, 
while Iqo and 7 remain the same. The tran- 
sient current therefore decays slower. This 
is because decreasing the couplings decreases 



the energy required for the particle to move 
around. This increase in the particle's free- 
dom to move increases the frequency of oscil- 
lation and slightly lengthens the decay of the 
transient current. However, fixing the value 
of V and increasing v^^ decreases ao and /qo, 
but increases /3. Therefore, for a given bias 
potential, increasing v^^ lengthens the decay 
but suppresses the amplitude of oscillation 
of the transient current. The parameters o^o, 

/qo, and 7 depend on the type of mate- 
rial used. The value of the interleads cou- 
pling v^^^ in addition, depends on the dis- 
tance between the leads. The farther apart 
are the two leads, the higher is the value of 
v^^ because of the higher interleads poten- 
tial barrier. Our results show that stronger 
interleads scattering lengthens the decay of 
the transient current. However, the scatter- 
ing also suppresses the amplitude of the tran- 
sient current and decreases the eventual value 
of the long-time steady-state current. 

The times when the peaks in the current 
occur can be known from the extremum of 
the power-law form of the current in Eq. f fTTl ) . 
Taking the time derivative of Eq. f fTTl ) and 
then equating the result to 0, we find the ex- 
tremum of the current to occur at times tp 
whenever the following is satisfied: 

= tan(cj + 0) . (14) 

a 

The left-hand side is an equation for a 
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TABLE II. (Color online) Values from the 
power-law fits to the exponent a as a function of 
the bias potential Ui, when the switch-on of the 
interleads coupling is in the form of a Heaviside 
step function. The shaded entries indicate cases 
when V ^ v^^. The dimension of is (1/eV)^ 
and /qo is /xA/(eV)^. 
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FIG. 6. (Color online) The location of the peaks 
in the transient current. The (black) lines repre- 
sent tan(a; t + 0), where uj — 2 [1/t] and = tt. 
The (red) dashed line represents uj t/a^ where 
a = 0.2. The (blue) dashed-dot line is when 
a = 0.3 and the (green) dashed-double dot line is 
when a — 0.5. The times when the straight and 
tangent lines intersect correspond to the times 
when the peaks in the transient current occur. 

straight line with a slope that depends on 
a. Since a varies depending on the values 



of the bias potential, the couphngs, and the 
speed of the switch-on, changing these pa- 
rameters would change the slope. As a con- 
sequence, the location in time of the current 
peaks would also change. This can be seen 
by noting how the peaks in the transient cur- 
rent move in Fig. [2] as a is varied, tp can 
be determined by the intersection points of 
the straight and tangent lines, corresponding 
to the left-hand side and right-hand side, re- 
spectively, of Eq. ([Hj) and as shown in Fig. [61 
The times when the current peaks occur are 
located whenever the two curves intersect. 
Since the slope of the straight line depends on 
a, we see that the faster decaying transient 
current, i.e., higher values of a, correspond 
to earlier peak times. 

Finally, we investigate the transport prop- 
erties of the device having a regular time 
variation, such as a nano-oscillator. In a 
nano-oscillator, the rotating disk in Fig. [l](a) 
is rocked back and forth across the dashed 
line. This would result in a harmonic modu- 
lation of the interleads coupling and would 
dynamically modulate the current through 
the device. However, compared to an al- 
ternating current which changes sign, the 
modulated current flowing out of the nano- 
oscillator maintains the same sign. We model 
the harmonically modulated coupling in the 
form v^^{t) = • (1 — coscj^it), where 

is the driving frequency of modulation. Fig- 
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FIG. 7. (Color online) The current as a func- 
tion of time for the nano-oscillator with driving 
frequency Ud — 0.25 [1/t]. (a) The bias poten- 
tials are ^7^ = 0.1 eV (black line), ^7^ = 0.3 eV 
(red line), and ^7^ = 0.5 eV (green line). The 
couplings are v = = — 2.7 eV. (b) The bias 
potential is = 0.5 eV and the hopping pa- 
rameter is V = —2.7 eV. The interleads coupling 
has amplitudes v^^ = —2.1 eV (black line) and 
v^^ = — 3.0 eV (orange line). The (blue) dashed 
lines show the harmonic variation of the inter- 
leads coupling. Their amplitudes are not drawn 
to scale. 

ure [71 shows the current characteristics as a 
function of time as the interleads coupling 
is swinged back and forth with driving fre- 
quency ujd = 0.25 [1/t]. In Fig. [7](a), the 
couplings are v = v^^ = —2.7 eV and the 
bias potential is varied. In Fig. [7](b) the bias 
potential is set at [/^ = 0.5 eV, the hopping 
parameter is flxed at = — 2.7 eV, and the 
interleads coupling v^^ is varied. We flnd 
that the current through the oscillator comes 



in pulses. However, it does not exactly follow 
the harmonic form of the coupling. The in- 
terleads coupling is maximum at times when 
the left and right leads are exactly aligned. 
In contrast, we see that the peaks in the cur- 
rent do not coincide with the times when the 
interleads coupling is maximum. The shape 
of the curve for the current actually looks like 
the truncated version of the transient current 
we examined in Fig. [2l In particular, the ini- 
tial overshoot of the transient current mani- 
fests as the first current peak in Fig. [71 This 
peak, however, does not occur when the in- 
terleads coupling is maximum. In addition, 
the times when the peak occurs depend on 
the values of the bias potential and the cou- 
plings. This dependence of the peak loca- 
tion to the above physical parameters follow 
similar dependence of the peak location in 
the nano-relay. In the design of nano-circuits 
containing an oscillator, therefore, it should 
be noted that the maximum current does not 
occur when the leads are exactly aligned and 
that the exact location of these peaks depend 
on the values of the applied bias potential, 
the nearest-neighbor coupling, and the inter- 
leads coupling. 



IV. SUMMARY AND 
CONCLUSION 

In summary, we examine a device that 
could act as a nano-relay or a nano-oscillator. 
The device consists of two leads and a time- 
varying interleads coupling. We use NEGF 
to derive a nonperturbative expression for the 
time-dependent current flowing from one lead 
to the other. In the nano-relay conflguration, 
we model the switch-on of the interleads cou- 
pling in the form of either a step function or 
a slowly progressing hyperbolic tangent. We 
flnd that the current oscillates and decays 
in time just after switch-on and during the 
transient regime. In both the step function 
and hyperbolic tangent switch-on, the de- 
cay of the transient current flts a power law. 
This leads to an equivalent RLC series circuit 
where all of the components have dynamical 
properties. We also flnd that the values of the 
couplings V and v^^^ the scattering at the in- 
terleads coupling, the speed of the switch-on, 
and the value of the bias potential affect the 
decay time of the transient current. In the 
long-time regime, the current approaches the 
steady-state value. In the nano-oscillator, we 
model the dynamical sytem by harmonically 
modulating the interleads coupling. We flnd 
that the current passes through the device in 
pulses, maintains the same sign, but does not 
exactly follow the functional form of the os- 



14 



cillating coupling. In particular, the peaks in 
the current do not occur at the times when- 
ever the leads are exactly aligned. 

The expressions for the current shown in 
Eq. (Ill) and the corresponding lesser Green's 
function shown in Eq. ([6j) are general and 
should be applicable to transport in quasi- 
linear systems where a switch-on in time oc- 
curs. The current oscillates and decays as 
a power law after a switch-on. This power- 
law decay implies the presence of dynamical 
resistance, inductance, and capacitance com- 
ponents. 
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